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' ABSTRACT 

O 

o ■ 

We present a sample of 17 high-redshift (3.5 ^ z ^ 5.2) QSO candidates in 
the 320 arcmin 2 area of the Great Observatories Origins Deep Survey, selected in 
C/j \ the magnitude range 22.45 < z 8 50 < 25.25 using deep imaging with the Advanced 

Camera for Surveys onboard the Hubble Space Telescope and the Advanced CCD 
! Imaging Spectrometer onboard the Chandra X-ray Observatory. On the basis 

of seven spectroscopic and ten photometric redshifts we estimate that the final 
sample will contain between two and four QSOs with 4 < z < 5.2. A dearth of 
high-redshift, moderate-luminosity (M 145 ~ —23) QSOs is observed with respect 
to predictions based on a) the extrapolation of the z ~ 2.7 luminosity function 
(LF), according to a pure luminosity evolution calibrated by the results of the 
\ Sloan Digital Sky Survey; and b) a constant universal efficiency in the formation 

Cp-T of super-massive black holes (SMBHs) in dark-matter halos. Evidence is gathered 

in favor of a density evolution of the LF at high redshift and of a suppression of 
' the formation or feeding of SMBHs in low-mass halos. 
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1. Introduction 

QSOs are intrinsically luminous and therefore can be seen rather easily at large dis- 
tances; however, they are rare, and so finding them requires surveys over large areas. As a 
consequence, at present, the number density of QSOs at high redshift is not well known. Re- 
cently, the Sloan Digital Sky Survey (SDSS) has produced a breakthrough, discovering QSOs 
up to z = 6.43 (Fan et al. 2003) and building a sample of six QSOs with z > 5.7. The SDSS, 
however, is sensitive only to very luminous QSOs (Mi 45 ^ —26.5) and provides no informa- 
tion about the faint end of the high- z QSO LF, which is particularly important to understand 
the interplay between the formation of galaxies and super-massive black holes (SMBH) and 
to measure the QSO contribution to the UV ionizing background (Madau, Haardt & Rees 
1999). New deep mult i- wavelength surveys like the Great Observatories Origins Deep Sur- 
vey (GOODS) provide significant constraints on the space density of less luminous QSOs 
at high redshift. Here we present high-z QSO candidates, identified in the two GOODS 
fields on the basis of deep imaging in the optical (with HST) and X-ray (with Chandra), 
and discuss the allowed space density of QSOs in the early universe. We adopt a definition 
of QSOs comprising all objects with strong, high-ionization emission lines and Mi 45 < — 21, 
including both conventional, broad-lined (type-1) QSOs and narrow-lined, obscured (type-2) 
QSOs. Throughout this paper we use a cosmology with h, f2 to t, — 0.7,1.0,0.3,0.7; 

magnitudes and colors are measured in the AB system. 

2. The Database 

The present work surveys an area of 320 arcmin 2 , subdivided in two fields centered on the 
Chandra Deep Field-South (CDF-S) and Hubble Deep Field-North (HDF-N), each covering 
m 10' x 16'. The optical data have been obtained with the Advanced Camera for Surveys 
(ACS) onboard HST in the framework of the GOODS/ACS survey described in Giavalisco 
et al. (2003a). Mosaics have been created from the first three epochs of observations, out 
of a total of five, in the bands F435W( J B 435 ), F606W(V r 606 ), F775W (i 775 ), F850LP(2 850 ). 
The catalogs used to select high-z QSOs have been generated using the SExtractor software, 
performing the detection in the z$ 50 band and then using the isophotes defined during this 
process as apertures for photometry in the other bands (Giavalisco et al. 2003a). 

The HDF-N and CDF-S fields have X-ray observations of 2 Ms and 1 Ms, respectively 



1 Based on observations taken with the NASA/ESA Hubble Space Telescope, which is operated by the 
Association of Universities for Research in Astronomy, Inc. (AURA) under NASA contract NAS5-26555. 
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(Alexander et al. 2003, hereafter A03; Giacconi et al. 2002, hereafter G02), providing 
the deepest views of the Universe in the 0.5-8.0 keV band. The X-ray completeness limits 
over ps 90% of the area of the GOODS fields are similar, with flux limits (S/N= 5) of 
~ 1.7 x 10~ 16 erg cm" 2 s" 1 (0.5-2.0 keV) and ps 1.2 x 10~ 15 erg cm" 2 s" 1 (2-8 keV) in the 
HDF-N field, and ps 2.2 x 10~ 16 erg cm" 2 s" 1 (0.5-2.0 keV) and ps 1.5 x 10~ 15 erg cm" 2 s" 1 
(2-8 keV) in the CDF-S field (A03). The sensitivity at the aim point is about 2 and 4 times 
better for the CDF-S and HDF-N, respectively. As an example, assuming an X-ray spectral 
slope of T = 2.0, a source detected with a flux of 1.0 x 10~ 16 erg cur 2 s _1 would have both 
observed and rest-frame luminosities of 8 x 10 42 erg s -1 , and 3 x 10 43 erg s -1 at z — 3, and 
z — 5, respectively (assuming no Galactic absorption). A03 produced point-source catalogs 
for the HDF-N and CDF-S and G02 for the CDF-S. 



3. Selection of the QSO candidates and determination of the redshifts 

The selection of the QSO candidates has been carried out in the magnitude interval 
22.45 < z 8 5o < 25.25. QSO colors in the ACS bands (Fig. 1) have been estimated as a 
function of redshift using a template (Cristiani & Vio 1990, CV90) of the QSO spectral 
energy distribution (SED) convolved with the Madau, Haardt & Rees (1999) model of the 
intergalactic medium (IGM) absorption. Four optical criteria have been tailored in order to 
select QSOs at progressively higher redshift in the interval 3.5 ^ z ^ 5.2: 

i-z < 0.35 AND 1.1 < B - V < 3.0 AND V - i < 1.0 (1) 

i-z < 0.35 AND B - V > 3.0 (2) 

i-z < 0.5 AND B - V > 2.0 AND V - i > 0.8 (3) 

i-z < 1.0 AND V - i > 1.9 (4) 

They select a broad range of high-z AGN, not limited to broad-lined (type-1) QSOs, and 
are less stringent than those typically used to identify high-z galaxies (e.g., Giavalisco et al. 
2003b). Below z ~ 3.5 the typical QSO colors in the ACS bands move close to the locus of 
stars and low-redshift galaxies. Beyond z ~ 5.2 the % — z color starts increasing and infrared 
bands would be needed to identify QSOs efficiently with an "i-dropout" technique. To avoid 
contamination from spurious sources we have limited our selection to zgso detections with 
S/N > 5. The criteria have been applied independently and produced in total 645 candidates 
in the CDF-S and 557 in the HDF-N. 

Two QSOs were found in the literature within the magnitude and redshift ranges of 
interest (22.45 < z 850 < 25.25, 3.5 < z < 5.2): CDF-S 033229.8 - 275105(z = 3.700) and 
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HDF-N 123647.9 + 620941(z = 5.186). Both have been selected with the present criteria 12 . 
Ten galaxies with 3.5 < z < 5.2 are known in the HDF-N (Cohen et al. 2000; Dawson et 
al. 2001, 2002). Eight are selected with the color criteria (1-4) (see Fig. 1). The two missed 
objects (HDF-N 1236279+6217504 and 1236376+621453) are both positionally uncertain 
identifications with a spectrum "solo-line with no continuum" (Dawson et al. 2001). HDF- 
N 1236279+6217504 has no optical counterpart in our ACS photometric catalog, and the 
supposed counterpart of HDF-N 1236376+621453 has colors very far from our criteria (-B435 — 
Vgog — 1-4, i 775 — z S50 = 1.3, 1^606 — ^775 = 2.3). 

The optical candidates selected with the criteria (1-4) have been matched with X-ray 
sources detected by Chandra (A03, G02) within an error radius corresponding to the 3 a 
X-ray positional uncertainty. With this tolerance the expected number of false matches 
is five and indeed two misidentifications, i.e. cases in which a brighter optical source lies 
closer to the X-ray position, have been rejected (both in the CDF-S). The sample has been 
reduced in this way to 11 objects in the CDF-S and 6 in the HDF-N (Tab. 1). Type-1 QSOs 
with M145 < —21, given the measured dispersion in their optical-to- X-ray flux ratio (Vignali 
et al. 2003), should be detectable in our X-ray observation up to z <^ 5.2. Conversely, 
any z > 3.5 source in the GOODS region detected in the X-rays must harbor an AGN 
(L x (0.5 - 2 keV) ^ 10 43 erg s" 1 ). 

Photometric redshifts of the 17 QSO candidates (Tab. 1, Col. 11) have been estimated 
by comparing with a x 2 technique (see Arnouts et al. (1999) for details) the observed ACS 
colors to those expected on the basis of a) a library of template SEDs of galaxies (the 
"extended Coleman" of Arnouts et al. (1999)); b) the QSO SED described in § 3. In general 
the estimates a) and b) are similar, since the criteria (1-4) both for galaxies and QSOs rely 
on a strong flux decrement in the blue part of the spectrum - due to the IGM and possibly 
an intrinsic Lyman limit absorption - superimposed on an otherwise blue continuum. Seven 
objects out of the 17 selected have spectroscopic confirmations (Tab. 1, Col. 12). Four are 
QSOs with redshifts between 2.7 and 5.2, in good agreement with the photometric redshifts 
based on the QSO SED. Three are reported to be galaxies and the relatively large offsets 
between the X-ray and optical positions suggest that they could be misidentifications. 



12 The X-ray emitting radio galaxy HDF-N 123642.0 + 621331, identified by Waddington et al. (1999) at 
z = 4.424, has not been selected. It should be noted, however, that the redshift of this source has been 
questioned by Barger, Cowie & Richards (2000). 
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Discussion and Conclusions 



We have compared the QSO counts observed in the zg5o band with two phenomeno- 
logical and two more physically motivated models (Fig. 2). In the absence of a complete 
spectroscopic follow-up we focus the comparison in the redshift range 4 < z < 5.2, where the 
selection criteria (1-4) and the photometric redshifts are expected to be highly complete and 
reliable. In addition to one spectroscopically confirmed QSO (HDF-N 123647.9 + 620941, 
z = 5.186) we estimate that between 1 — 3 more QSOs at z > 4 are present in the GOODS, 
depending on whether galaxy or QSO SEDs are adopted for the photometric redshifts. 

Phenomenological models. The double power-law fit of the 2QZ QSO LF (Boyle et al. 
2000) has been extrapolated for z > 2.7 (the peak of QSO activity) in a way to produce a 
power-law decrease of the number of bright QSOs by a factor 3.5 per unit redshift interval, 
consistent with the 3.0*g g factor found by Fan et al. (2001). The extrapolation is carried out 
either as a Pure Luminosity Evolution (PLE) or a Pure Density Evolution (PDE), assuming 
that the slopes of the LF power-laws remain unchanged all the way up to z ~ 5. The PLE 
model predicts about 17 QSOs with z 850 < 25.25 at redshift z > 4 (27 at z > 3.5) in the 320 
arcmin 2 of the two GOODS fields, and is inconsistent with the observations at a more than 
3a level. The PDE estimate is 2.9 QSOs at z > 4 (6.7 at z > 3.5). 

Physically motivated models. It is possible to connect QSOs with DMH formed in 
hierarchical cosmologies with a minimal set of assumptions (MIN model, e.g. Haiman & Hui 
(2001)): a) QSOs are hosted in newly formed halos with b) a constant SMBH/DMH mass 
ratio e and c) accretion at the Eddington rate. The bolometric LF of QSOs is then: 



where the abundance of DMHs of mass M H is computed using the Press & Schechter (1974) 
recipe, the distribution P(tf\M H) t(z))dtf of the formation times tf follows Lacey & Cole 
(1993), tduty is the QSO duty cycle and Lem = 10 4 53 L & is the Eddington luminosity of a 
1 M Q SMBH. The MIN model is known to overproduce the number of high- z QSOs (Haiman, 
Madau & Loeb 1999) and in the present case predicts 151 QSOs at z > 4 (189 at z > 3.5). 
Alternatively, QSOs are assumed to shine a tdciay time after DMH formation (DEL model; 
Monaco, Salucci & Danese 2000; Granato et al. 2001). The QSO LF of eq. 5 is then 
computed at a redshift z' corresponding to t(z) — tdciay The predictions of the DEL model 
are very close to the PDE, with 3.2 QSOs expected at z > 4 (4.8 at z > 3.5), in agreement 
with the observations. While all the models presented here are consistent with the recent 
QSO LF measurement of the COMBO-17 survey (Wolf et al. 2003) at z ~ 2 - 3, only the 
PDE and DEL models fit the COMBO-17 LF in the range 4.2 < z < 4.8 and M U5 < -26. 
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At z > 4 the space density of moderate luminosity (M U5 ~ —23) QSOs is significantly 
lower than the prediction of simple recipes matched to the SDSS data, such as a PLE 
evolution of the LF or a constant universal efficiency in the formation of SMBH in DMH. 
A similar result has been obtained at 5 ^ z ^ 6.5, by Barger et al. (2003) who also showed 
that the QSO contribution to the UV background is insufficient to ionize the IGM at these 
redshifts. A flattening of the observed high- z LF is required below the typical luminosity 
regime (M 145 ^ —26.5) probed by the SDSS. This is an indication that at these early epochs 
the formation or the feeding of SMBH is strongly suppressed in relatively low- mass DMH, as 
a consequence of feedback from star formation (Granato et al. 2001) and/or photoionization 
heating of the gas by the UV background (Haiman, Madau & Loeb 1999). 

We are grateful to H.Spinrad, S.Dawson, and A.Dey for providing the redshift for HDF-N 
123627.5 + 621158 and to C.Steidel for the spectrum of the galaxy HDF-N 123644.1 + 621311. 
The work of DS was carried out at the Jet Propulsion Laboratory, California Institute of 
Technology, under a contract with NASA. DMA, FEB, and WNB thank the NSF CAREER 
award AST-9983783. DMA also acknowledges the support provided by a Royal Society 
University Research Fellowship. We thank STScI grant HST-GO-09425.26-A (FEB, WNB), 
HST-GO-09425.13-A (CMU, EC, ET), and ASI grant I/R/088/02 (SC, AG, MN, EV). 
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Table 1. High-redshift QSO candidates in the CDF-S. 
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'Objects with upper limits both in the soft and hard band have been detected in other sub-bands. 
§The first photometric redshift is computed on the basis of galaxy SEDs, the second for QSO SEDs. 
+ From the supplementary X-ray catalog by A03. 

♦Detected in the X-ray by G02 and not by A03. X-ray coordinates and fluxes derived from G02. Positions in G02 arc typically 
offset +0.07" in RA and —1" in DEC with respect to the optical ACS positions. The correction has been applied in Cols. 3-4. 
"QSO showing relatively narrow emission lines with a P-Cyg profile (Cristiani et al. 2000). 
b Type-2 QSO (Norman ct al. 2002). 
°Galaxy: H. Spinrad, private communication. 
d Broad emission-line QSO (Bargcr ct al. 2001). 
c Galaxy identified by Cohen et al. (2000). 

f Galaxy identified by Dickinson (1998). No spectroscopic evidence for AGN activity is present in the spectrum. 
«QSO identified by Barger et al. (2002). 
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Fig. 1. — Expected 2775 — zsso vs. -B435 — Veoe colors of QSOs as a function of redshift 
compared to other classes of objects. The dashed line shows the locus of QSOs at z > 3.5 
estimated with the composite spectrum of CV90. The continuous and dotted lines show 
the colors of Ly-break and Sy 2 galaxies, derived from the SEDs of Schmitt et al. (1997) 
and Arnouts et al. (1999), respectively. The corresponding symbols - open triangles, circles 
and squares - start on the left at z — 3.5 and move to the right in steps of Az = 0.1. The 
observed colors of five QSOs are marked with filled circles and the corresponding redshifts. 
The same is done for ten high- ,2 galaxies identified with crosses. The radio galaxy HDF-N 
123642.0 + 621331 is marked with the label 4.4W. The positions of various types of stars is 
also shown. The dot-dashed line represents the projection of the selection criteria 1) + 3). 
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Fig. 2. — Comparison of the observed QSO counts and redshift distribution with model 
predictions. Circles and stars show the "maximal" and "minimal" estimates of the GOODS 
counts, respectively (see text). The dotted segments show the corresponding 1 a upper 
(maximal case) and lower (minimal case) confidence limits. Models are represented by four 
smooth lines: dot-dashed PLE, continuous PDE, short-dash MIN (Haiman & Hui 2001), 
long-dash delayed QSO shining. A value of a% = 0.83 has been assumed in the MIN and 
DEL models, whose parameters have been fixed in order to reproduce the QSO LF at z ~ 3. 
In the MIN model: loge = -3.2 and t duty = 0.015 Gyr. In the DEL model: loge = -3 
and tduty = 0.013 Gyr; the delay time tdeiay(Mff) is set to 0.75 Gyr above 10 125 M® and 
is decreased as cx M^ ' 7 at lower masses (Granato et al. 2001). As dynamical and cooling 
times are shorter at high redshift, tdeiay is assumed to decrease as (1 + z)~ 2 above z = 2.7. 



